Laser cooling of a nanomechanical oscillator into its quantum ground state 
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A patterned Si nanobeam is formed which supports co-localized acoustic and optical resonances that are 
coupled via radiation pressure. Starting from a bath temperature of Tfy ^ 20 K, the 3.68 GHz nanomechanical 
mode is cooled into its quantum mechanical ground state utilizing optical radiation pressure. The mechanical 
mode displacement fluctuations, imprinted on the transmitted cooling laser beam, indicate that a final phonon 
mode occupancy ofn = 0.85 ± 0.04 is obtained. 



The simple mechanical oscillator, canonically consisting of a coupled mass-spring system, is used in a wide variety of sensitive 
measurements, including the detection of weak forces 1 1 1 and small masses |2|. A classical oscillator can take on a well-defined 
amplitude of sinusoidal motion. A quantum oscillator, however, has a lowest energy state, or ground state, with a finite amplitude 
uncertainty corresponding to the zero-point motion. In our everyday experience mechanical oscillators are filled with many 
energy quanta due to interactions with their highly fluctuating thermal environment, and the oscillator's quantum nature is all 
but hidden. Recently, in experiments performed at temperatures of a few hundredths of a Kelvin, engineered nanomechanical 
resonators coupled to electrical circuits have been measured to be oscillating quietly in their quantum ground state |T, 4|. These 
experiments, in addition to providing a glimpse into the underlying quantum behavior of mesoscopic systems consisting of 
billions of atoms, represent the initial steps towards the use of mechanical elements as tools for quantum metrology |[5]|6l or as 
a means to couple hybrid quantum systems fT-^ . In this work we have created a coupled, nanoscale optical and mechanical 
resonator 1 10] formed in a silicon microchip, in which radiation pressure from a laser is used to cool the mechanical motion 
down to the quantum ground state (average phonon occupancy, /2 = 0.85 ±0.04). Critically, this cooling is realized at an 
environmental temperature some thousand times larger than in previous experiments (7], ^ 20 K), and paves the way for optical 
control of mesoscale mechanical oscillators in the quantum regime. 

It has been known for some time 1 1 1 1 that atoms and ions nearly resonant with an applied laser beam (or series of beams) may 
be mechanically manipulated, even trapped and cooled down to the quantum ground state of their center-of-mass motion 1 12 |. 
Equally well known (Tl has been the fact that radiation pressure can be exerted on regular dielectric (i.e., non-resonant) objects 
to damp and cool their mechanical motion. In so-called cavity-assisted schemes, the radiation pressure force is enhanced by 
coupling the motion of a mechanical object to the light field in an optical cavity. Pumping of the optical cavity by a single- 
frequency electromagnetic source produces a coupling between the mechanical motion and the intensity of the electromagnetic 
field built-up in the resonator. As the radiation pressure force exerted on the mechanical object is proportional to the field 
intensity in the resonator, a form of dynamical back-action results dllSl- For a lower- frequency (red) detuning of the laser from 
the cavity, this leads to damping and cooling of the mechanical motion. 

Recent experiments involving micro- and nanomechanical resonators, coupled to electromagnetic fields at optical and mi- 
crowave frequencies, have demonstrated significant radiation pressure dynamic back-action 1 13 1. These structures have included 
Fabry-Perot cavities with mechanically-compliant miniature end-mirrors [il4UT6l or internal nanomembranes IfTi . whispering- 
gallery glass resonators |18|, nanowires capacitively coupled to co-planar microwave transmission-line cavities |[6j[l9|, and 
lumped circuit microwave resonators with deformable, nanoscale, vacuum-gap capacitors 1201 . The first measurement of an 
engineered mesoscopic mechanical resonator predominantly in its quantum ground state, however, has been performed not us- 
ing back-action cooling, but rather, using conventional cryogenic cooling (bath temperature ^ 25 mK) of a high frequency, 
and thus lower thermal occupancy, oscillator O. Read-out and control of mechanical motion at the single quanta level was 
performed by strongly coupling the GHz-frequency piezoelectric mechanical resonator to a resonant superconducting quantum 
circuit. Only recently have microwave systems, also operating at bath temperatures of 7], ^ 25 mK, utilized radiation pressure 
back-action to cool a high-g, MHz-frequency mechanical oscillator to the ground state | 4, 19 1. 

Optically coupled mechanical devices, while allowing for control of the mechanical system through well-established quantum 
optical techniques |21 1, have thus far not reached the quantum regime due to a myriad of technical difficulties 1 18 |. A particular 
challenge has been maintaining efficient optical coupling and low-loss optics and mechanics in a cryogenic, sub-Kelvin environ- 
ment. The optomechanical system studied in this work enables large optical coupling to a high- 2 GHz-frequency mechanical 
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FIG. 1: Optomechanical resonator with phononic shield, a, Scanning electron microscope (SEM) image of the patterned Si nanobeam 
with external phononic bandgap shield, b, Enlarged SEM image of the central cavity region of the nanobeam. c, Top, finite element method 
(FEM) simulated normalized electric field of the localized optical resonance of the nanobeam cavity. Bottom, FEM simulated normalized 
displacement field of the acoustic resonance (breathing mode) which is coupled via radiation pressure to the co-localized optical resonance. 
The displacement field is indicated by the exaggerated deformation of the structure, with the relative magnitude of the local displacement 
(strain) indicated by the color, d, SEM image of the interface between the nanobeam and the phononic bandgap shield, e, FEM simulation 
of the normalized squared displacement field amplitude of the localized acoustic resonance at the nanobeam- shield interface, indicating the 
strong suppression of acoustic radiation provided by the phononic bandgap shield. The color scale represents log(x^/ max{x^}) where x is the 
displacement. 



oscillator, allowing for both efficient back-action cooling and significantly higher operating temperatures. As shown in Fig.[T^, 
the system consists of an integrated optical and nanomechanical resonator formed in the surface layer of a silicon-on-insulator 
microchip. The periodic patterning of the nanobeam is designed to result in Bragg scattering of both optical and acoustic guided 
waves. A perturbation in the periodicity at the center of the beam results in co-localized optical and mechanical resonances 
(Fig. [T]3-c), which are coupled via radiation pressure ifTOl . The fundamental optical resonance of the structure occurs at a fre- 
quency C0o/27r = 195 THz (k = 1537 nm), while, due to the much slower speed of sound, the mechanical resonance occurs at 
co^/27r = 3.68 GHz. In order to minimize mechanical damping in the structure, an external acoustic radiation shield is added in 
the periphery of the nanobeam (Fig.[T]l-e). This acoustic shield consists of a two-dimensional "cross" pattern, which has been 
shown both theoretically and experimentally to yield a substantial phononic bandgap in the GHz frequency band f22l. 

A fiber taper nanoprobe, formed from standard single-mode optical fiber, is used to optically couple to the silicon nanoscale 
resonators. As shown in Fig.[2j a tunable laser (New Focus Velocity swept laser; 200 kHz linewidth) is used to optically cool and 
transduce the mechanical motion of the nanomechanical oscillator. Placing the optomechanical devices into a continuous-flow 
helium cryostat provides a modicum of pre-cooling down to Tij ^ 20 K, reducing the bath occupancy of the 3.68 GHz mechanical 
mode to nij ^ 100. At this temperature the mechanical Q-f actor increases up to a measured value of 2^ ^ 10^, corresponding to 
an intrinsic mechanical damping rate of yt/ln = 35 kHz. The optical Q-factor is measured to be 2^ = 4 x 10^, corresponding to 
an optical linewidth of K/27r = 500 MHz, slightly reduced from its room temperature value. 

In the resolved sideband limit in which co^/k > 1, driving the system with a laser (frequency CO/) tuned to the red side of the 
optical cavity (detuning A = co^ — CO/ = co^), creates an optically-induced damping, yoM, of the mechanical resonance |23, 24l . 
In the weak-coupling regime (yoM ^ k) and for a detuning A = co^, the optical back-action damping is given by yoM = 4^^^c/k, 
where ric is the average number of drive photons stored in the cavity and g is the optomechanical coupling rate between the 
mechanical and optical modes. This coupling rate, is quantified as the shift in the optical resonance for an amplitude of 
motion equal to the zero-point amplitude (Xzpf = (fi/2mco^)^/^; m the motional mass of the localized acoustic mode, h Planck's 
constant divided by In). The optomechanical damping, a result of the preferential scattering of drive laser photons into the upper- 
frequency sideband, also cools the mechanical mode. For a quantum-limited drive laser, the phonon occupancy of the mechanical 
oscillator can be reduced from nij = ksTij/fKi^m ^ 1, to a value /r = ^/^/(l + C) + ^min, where C = Jou/ji is the cooperativity. 
The residual scattering of drive laser photons into the lower-frequency sideband limits the cooled phonon occupancy to ^min = 
(k/4co^)^, determined by the level of sideband resolution 12311241 . 

The drive laser, in addition to providing mechanical damping and cooling, can be used to measure the mechanical and optical 
properties of the system through a series of calibrated measurements. In a first set of measurements the noise power spectral 
density (PSD) of the drive laser transmitted through the optomechanical cavity is used to perform spectroscopy of the mechanical 




FIG. 2: Experimental setup. A single tunable 1550 nm diode laser is used as the cooling and mechanical transduction beam sent into the 
nanobeam optomechanical resonator cavity held in a continuous flow Helium cryostat. A wavemeter (WM) is used to track and lock the laser 
frequency, while a variable optical attenuator (VOA) is used to set the laser power. The transmitted signal is amplified by an erbium doped fiber 
amplifier (EDFA), and detected on a high-speed photodetector (D2) connected to a real-time spectrum analyzer (RSA), where the mechanical 
noise power spectrum is measured. A slowly modulated probe signal used for optical spectroscopy and calibration is generated from the 
cooling laser beam via an amplitude electro-optic modulator (EOM) driven by a microwave source (RFSG). The reflected component of this 
signal is separated from the input via an optical circulator, sent to a photodetector (Dl), and then demodulated on a lock-in amplifier. Paddle- 
wheel fiber polarization controllers (FPCs) are used to set the laser polarization at the input to the EOM and the input to the optomechanical 
cavity. For more detail see Appendix [P] 



mode. As shown in Appendix |A| the noise power spectral density of the photocurrent generated by the transmitted field of 
the drive laser with red-sideband detuning (A = co^) yields a Lorentzian component of the single-sided PSD proportional to 
Sh (co) = ny/ ( (co — co^)^ + (y/2)^) , where 7 = 7;+ 7om = 7- ( 1 + C) is the total mechanical damping rate. For a blue laser detuning 
of A = — cOm, the optically-induced damping is negative (70M = — 4g^^c/K) and the photocurrent noise PSD is proportional 
to 5^1(03) = {n-\- 1)7/ ((co — co^)^ + (7/2)^). Typical measured noise power spectra under low power laser drive (ric = 1.4, 
C = 0.27), for both red and blue detuning, are shown in Fig.|3^. Even at these small drive powers the effects of back-action are 
clearly evident on the measured spectra, with the red-detuned drive broadening the mechanical line and the blue-detuned drive 
narrowing the line. The noise floor in Fig. [3^ (shaded in gray) corresponds to the noise generated by the erbium doped fiber 
amplifier (EDFA) used to pre-amplify the transmitted drive laser signal prior to photodetection, and is many orders of magnitude 
above the electronic noise of the photoreceiver and real-time spectrum analyzer. 

Calibration of the EDFA gain, along with the photoreceiver and real-time spectrum analyzer photodetection gain, allows one 
to convert the measured area under the photocurrent noise PSD into a mechanical mode phonon occupancy. As described in 
detail in the Appendices, these calibrations, along with measurements of low drive power (C <C 1) rf-spectra of both A = ±co^ 
detunings, are performed to provide an accurate, local thermometry of the optomechanical cavity. An example of this form of 
calibrated mode thermometry is shown in Eig.jSj), where we plot the optically measured mechanical mode bath temperature (T^) 
versus the cryostat sample mount temperature (7^; independently measured using a Si diode thermometer attached to the copper 
sample mount). As one can see from this plot, the optical mode thermometry accurately predicts the absolute temperature of the 
sample for Tc > 50 K; below this value the mode temperature deviates from Tc and saturates to a value of 7], = 17.6±0.8K due 
to black-body heating of the device through the imaging aperture in the radiation shield of our cryostat. 

In a second set of measurements, the mechanical damping, 7, and the cavity-laser detuning. A, can be measured by optical 
spectroscopy of the driven cavity. By sweeping a second probe beam of frequency 0)^ over the cavity, with the cooling beam 
tuned to A = co^, spectra exhibiting electromagnetically-induced transparency (EIT) ||20l |25] |26l are measured, as shown in 
Fig. [3]:. Due to the high single-photon cooperativity in the system, an intracavity population of only ric ^ 5 switches the system 
from reflecting to transmitting for the probe beam. The corresponding dip at the center of the optical cavity resonance occurs 
at a two-photon detuning A^/ = co^ — CO/ = co^ and has a bandwidth equal to the mechanical damping rate, 7 (1 + C). Figure 
shows a plot of the measured mechanical linewidth versus intracavity photon number, displaying good correspondence between 



4 




Ay27i (GHz) 

FIG. 3: Mechanical and optical response, a, Typical measured mechanical noise spectra around the resonance frequency of the breathing 
mode for low laser drive power (ric = 1.4). The blue (red) curve corresponds to the measured spectrum with laser drive blue (red) detuned 
by a mechanical frequency from the optical cavity resonance. The black trace corresponds to the measured noise floor (dominated by EDFA 
noise) with the drive laser detuned far from cavity resonance, b, Plot of the measured (□) mechanical mode bath temperature (Tfy) versus 
cryostat sample mount temperature (Tc). The dashed line indicates the curve corresponding to perfect following of the mode temperature with 
the cryostat temperature (Ti^ = Tc). c, Typical reflection spectrum of the cavity while driven by the cooling laser at A = co^ as measured by a 
weaker probe beam at two-photon detuning A^/. The signature reflection dip on-resonance with the bare cavity mode, highlighted in the inset, 
is indicative of electromagnetically-induced transparency (BIT) caused by the coupling of the optical and mechanical degrees of freedom by 
the cooling laser beam. 



both mechanical and optical spectroscopy techniques. From a fit to the measured mechanical damping rate versus ric (dashed 
red line in Fig. [4^), the zero-point-motion optomechanical coupling rate is estimated to be g/27l = 910 kHz, placing the system 
well within the weak-coupling regime for all measured drive powers. 

In Fig.|4]3 we plot the calibrated Lorentzian noise PSD area, in units of phonon occupancy, versus red-detuned (A = co^) drive 
laser power. Due to the low effective temperature of the laser drive, the mechanical mode is not only damped but also cooled 
substantially. The minimum measured mode occupancy for the highest drive power of ric ^ 2000 is n = 0.85 ±0.04, putting 
the mechanical oscillator in a thermal state with ground state occupancy probability greater than 50%. The dashed blue line in 
Fig.Hb represents the ideal back-action cooled phonon occupancy estimated using both the measured mechanical damping rate in 
Fig.|^ and the low drive power intrinsic mechanical damping rate. Deviation of the measured phonon occupancy from the ideal 
cooling model is seen to occur at the highest drive powers, and as detailed in Appendices |I][K| is due to both an increase in the 
bath temperature due to optical absorption (Fig. [4]:) and an increase in the intrinsic mechanical damping rate (Fig.|4]l) induced 
by the generation of free-carriers through optical absorption. In order to evaluate the efficiency of the optical transduction of the 
mechanical motion, we also plot in Fig.|4^ the measured background noise PSD (or imprecision level) alongside that for an ideal 
cavity transducer with shot-noise limited detection. The minimum measurement imprecision corresponds to ^imp ^ 20 in units 
of phonon quanta, 40% of which is due to light lost inside the optical cavity (and thus not detected). The remaining 60% of the 
imprecision stems from optical loss in the fiber taper (approximately 2 dB) and added noise due to the EDFA pre-amplification. 

Looking ahead, the combination of strong optical back-action cooling and efficient mechanical motion transduction realized 
in the chip-scale optomechanical cavities of this work, represents a first step towards optical quantum control of nanomechanical 
objects |27|. For example, optomechanical entanglement between light and mechanics | 28 | or quantum state transfer between 
single optical photons and mechanical phonons 121 |29l may be envisioned, enabling mechanical systems to function as either 
quantum transducers fS] or quantum memory elements f30]. The efficacy of such quantum protocols and devices, and the ability 
to measure quantum dynamics, relies on the thermal decoherence time of the mechanical system, given by x = ksTiy/hQm- For 
the measured devices in this work at Tjy ^ 20 K, the decoherence time corresponds to A^osc = xco^/27r ^ 200 periods of coherent 
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FIG. 4: Optical cooling results, a, Measured mechanical mode linewidth (□), BIT transparency bandwidth (o), and predicted optomechanical 
damping rate estimated using the zero-point optomechanical coupling rate g/2n = 9\0 kHz (red dashed line). The inset shows the measured 
EIT transparency window at the highest cooling drive power, b, Measured (o) average phonon number, n, in the breathing mechanical mode 
at (Om/2n = 3.68 GHz versus cooling laser drive power (in units of intracavity photons, ric), as deduced from the calibrated area under the 
Lorentzian lineshape of the mechanical noise power spectrum, Si^. The left and right inset spectra correspond to the measured noise power 
spectrum in units of m^/Hz at low and high laser cooling power, respectively. The dashed blue line indicates the estimated mode phonon 
number from the measured optical damping alone. Error bars indicate computed standard deviations as outlined in the SI. c. Estimated bath 
temperature, T^, versus cooling laser intracavity photon number, ric. d. Measured change in the intrinsic mechanical damping rate versus ric 
(o). A polynomial fit to the mechanical damping dependence on ric is shown as a dashed line. For more detail see the SI. e. The measured 
(□) background noise PSD versus laser drive power (ric), in units of phonon quanta. Shown as a black dashed line is the theoretical shot-noise 
limited noise PSD for an ideal single-sided cavity, a unit quantum efficiency detector, and no optical loss in the transmitted optical field. 



oscillation of the mechanical resonator. Back-action cooling ll24ll from a given bath temperature to the quantum ground state 
requires Hk ^ ksTh/Qm. Given the properties of the devices in this work, ground state cooling from room temperature with 
measurable quantum dynamics (A/qsc ^10) seems feasible. This opens up the possibility for future experimentation with, and 
utilization of, quantum nanomechanical objects in a room temperature environment. 

Acknowledgements 

This work was supported by the DARPA/MTO ORCHID program through a grant from AFOSR, the European Commission 
(MINOS, QUESSENCE), European Research Council (ERC QOM), the Austrian Science Fund (CoQuS, FOQUS, START), and 
the Kavli Nanoscience Institute at Caltech. JC and ASN gratefully acknowledge support from NSERC. 



[1] Braginsky, V. & Manukin, A. Measurement of weak forces in Physics experiments (Univ. of Chicago Press, 1977). 

[2] Jensen, K., Kim, K. & Zettl, A. An atomic-resolution nanomechanical mass sensor. Nature Nanotech. 3, 533-537 (2008). 

[3] O'Connell, A. D. al Quantum ground state and single-phonon control of a mechanical resonator. Nature 464, 697-703 (2010). 

[4] Teufel, J. D. et al. Sideband cooling micromechanical motion to the quantum ground state. arXiv: 1103.2 144 (201 1). 

[5] Caves, C, Thome, K., Drever, R., Sandberg, V. D. & Zimmermann, M. On the measurement of a weak classical force coupled to a 

quantum-mechanical oscillator. Reviews of Modern Physics 52, 341-392 (1980). 
[6] Regal, C. A., Teufel, J. D. & Lehnert, K. W. Measuring nanomechanical motion with a microwave cavity interferometer. Nature Phys. 4, 

555-560 (2008). 

[7] Wallquist, M., Hammerer, K., Rabl, P., Lukin, M. & Zoller, P. Hybrid quantum devices and quantum engineering. Physica Scripta 2009, 
014001 (2009). 



6 



[8] Stannigel, K., Rabl, R, S0rensen, A. S., Zoller, R & Lukin, M. D. Optomechanical Transducers for Long-Distance Quantum Communi- 
cation. Phys. Rev. Lett. 105, 220501 (2010). 

[9] Safavi-Naeini, A. H. & Painter, O. Proposal for an optomechanical traveling wave phonon-photon translator. New J. Phys. 13, 013017 
(2011). 

[10] Eichenfield, M., Chan, J., Camacho, R. M., Vahala, K. J. & Painter, O. Optomechanical crystals. Nature 462, 78-82 (2009). 
[11] Cohen-Tannoudji, C. N. & Phillips, W. D. New Mechanisms for Laser Cooling. Phys. Today 43, 33-40 (1990). 

[12] Monroe, C. et al. Resolved-Sideband Raman Cooling of a Bound Atom to the 3D Zero-Point Energy. Phys. Rev. Lett. 75, 4011^014 
(1995). 

[13] Kippenberg, T. J. & Vahala, K. J. Cavity Optomechanics: Back- Action at the Mesoscale. Science 321, 1172-1176 (2008). 
[14] Gigan, S. et al. Self-cooling of a micromirror by radiation pressure. Nature 444, 67-70 (2006). 

[15] Arcizet, O., Cohadon, P.-R, Briant, T., Pinard, M. & Heidmann, A. Radiation-pressure cooling and optomechanical instabihty of a 

micromirror. Nature 444, 71-74 (2006). 
[16] Groblacher, S. et al. Demonstration of an ultracold micro-optomechanical oscillator in a cryogenic cavity. Nature Phys. 5, 485^88 

(2009) . 

[17] Thompson, J. D. et al. Strong dispersive coupling of a high-finesse cavity to a micromechanical membrane. Nature 452, 72-75 (2008). 
[18] Riviere, R. et al. Optomechanical sideband cooling of a micromechanical oscillator close to the quantum ground state. arXiv:101L0290 

(2010) . 

[19] Rocheleau, T. et al. Preparation and detection of a mechanical resonator near the ground state of motion. Nature 463, 72-75 (2010). 
[20] Teufel, J. D. et al. Circuit cavity electromechanics in the strong-coupling regime. Nature 471, 204-208 (2011). 
[21] Wiseman, H. M. & Milburn, G. J. Quantum Measurement and Control (Cambridge University Press, 2010). 

[22] Alegre, T. P. M., Safavi-Naeini, A., Winger, M. & Painter, O. Quasi-two-dimensional optomechanical crystals with a complete phononic 

bandgap. Opt. Express 19, 5658-5669 (2011). 
[23] Wilson-Rae, L, Nooshi, N., Zwerger, W. & Kippenberg, T. J. Theory of ground state cooling of a mechanical oscillator using dynamical 

back-action. Phys. Rev. Lett. 99, 093901 (2007). 
[24] Marquardt, R, Chen, J. P., Clerk, A. A. & Girvin, S. M. Quantum theory of cavity-assisted sideband cooling of mechanical motion. Phys. 

Rev. Lett. 99, 093902 (2007). 
[25] Weis, S. et al. Optomechanically Induced Transparency. Science 330, 1520-1523 (2010). 

[26] Safavi-Naeini, A. H. et al. Electromagnetically induced transparency and slow light with optomechanics. Nature 472, 69-73 (2011). 
[27] Aspelmeyer, M., Groblacher, S., Hammerer, K. & Kiesel, N. Quantum optomechanics - throwing a glance. / Opt. Soc. Am. B 27, 
A189-A197 (2010). 

[28] Vitah, D. et al. Optomechanical Entanglement between a Movable Mirror and a Cavity Pield. Phys. Rev. Lett. 98, 030405 (2007). 
[29] Akram, U., Kiesel, N., Aspelmeyer, M. & Milburn, G. J. Single-photon optomechanics in the strong coupling regime. New J. Phys. 12, 
083030 (2010). 

[30] Chang, D., Safavi-Naeini, A. H., Hafezi, M. & Painter, O. Slowing and stopping light using an optomechanical crystal array. New J. 
Phys. 13, 023003 (2011). 

[31] Borselli, M., Johnson, T. J. & Painter, O. Measuring the role of surface chemistry in silicon microphotonics. App. Phys. Lett. 88, 131114 
(2006). 

[32] Prey, B. J., Leviton, D. B. & Madison, T. J. Temperature-dependent refractive index of silicon and germanium. In Proc. SPIE, vol. 6273, 
62732J (2006). 

[33] Harrington, R. F. Time-Harmonic Electromagnetic Fields (IEEE Press, 1961). 

[34] Barclay, P., Srinivasan, K. & Painter, O. Nonlinear response of silicon photonic crystal microresonators excited via an integrated waveg- 
uide and fiber taper. Opt. Express 13, 801-820 (2005). 

[35] Desurvire, E., Bayart, D., Desthieux, B. & Bigo, S. Erbium-Doped Fiber Amplifiers, Device and System Developments (Wiley- 
Interscience, 2002). 



Appendix A: Classical Derivation of Transduced Signal 

We begin by modeling the optomechanical system with the Hamiltonian 

H = hAa^d-\-h(Ornb^b-\-hg{P +b)d^d-\-ihJ -^Oiin,o{^- d'^), (Al) 

where A = co^, — coi, with laser frequency coi, optical mode frequency (x)o and mechanical mode frequency co^. Here d (a^) and 
b (b^) are respectively the annihilation (creation) operators of photon and phonon resonator quanta, g is the optomechanical 
coupling rate corresponding physically to the shift in the optical mode frequency due to the zero-point fluctuations (xzpf = 
^yh/2m(0m, m motional mass) of the phonon mode. By making the substitutions 

d^a = Y, OLqe-'^"^^' , b po^"''*'"' (A2) 

we can treat the system classically by representing the photon amplitudes as a Fourier decomposition of sidebands. Notice that 
the infinite summation over each sideband order q, can be relaxed to a few orders in the sideband resolved regime (k <^ co^). The 
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phonon amplitude, po. is the classical mechanical excitation amplitude. For an oscillator undergoing thermal Brownian motion, 
po, is a stochastic process. We assert the stochastic nature of the variable, at the end of the derivation where the power spectral 
density is calculated. The equation of motion for the slowly varying component is then 



(A3) 



where we introduce the cavity (optical) energy loss rate, k, and the cavity coupling rate, k^. This can be written as a system of 
equations M • a = ain where 



M 



pq 



"■m,p 



-OCin,o8pO- 



(A4) 
(A5) 



By truncating and inverting the coupling matrix M one can determine each one of the sidebands amplitude as = (M ^)qp ain,p 
and therefore determine the steady state power leaving the cavity to be 



OCout — OCin,0 + a / 



(A6) 



where we assumed that the input pump is not depleted by the cavity, which is the frame of interest of this work. In this case the 
total power measured at the photodetector will be proportional to 



OCout = 



OCin,0 + A/ yOC 



q p [ ^ q J 



(A7) 
(A8) 



1. Resolved Sideband Limit 



The equations presented in the previous section are exact and therefore can be solved for any case. However our interests 
lie in the so-called resolved sideband limit, 0)^ > k/2, where further simplification can be done. Specifically for our system, 
co^/27r = 3.68 GHz, K/2n = 500 MHz putting us well within this limiting case. 

Additionally, in the cavities studied, the optomechanical phase-modulation factor (proportional to gXzpf/co^), is much less 
than 10~^. As such, only the 0) = COi zb co^ sidebands (^ = ±1) are significant and ao ^ a±. Truncating the matrix equations 
appropriately, we find 



oco = 



-v/K^ain,Q 
/A + k/2 

-^'gpooco 
/(A=fco^) + k/2 



(A9) 
(AlO) 



where ain,o = >/Mn, = P'm/^^o and Pin the input power at the cavity. 



In the sideband resolved limit, we have |ain,o| > | v^Kg/2ao| > |i/Kg/2a±|. Therefore the photodetector signal is predomi- 
nantly composed of the mixing between sidebands with the input pump beam and terms containing |aop, and can be written 
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as 



cni: 



icxoutp = |ain,oP + \l Yain,o(ao + aS) + y |ao|^ + ... 
2 



ain,o (a+ ^^^-^ + a; ^-^^-0 + ( I oco 1 1 a± I ) 



(All) 



|OCin,o| 



K./2 



/A + k/2 

cos(co^r) |A+|cos((p+) + |A_|cos((p_) 
sin(co^^) |A+|sin((p+)-|A_|sin((p_) 



(A12) 



where A± = 2^/K^/2 ain,oOC± = |A± | exp (— /(p±). We can easily recognize the first term in Equation ( A12 ) as the DC cavity trans- 
mission spectra. The remaining two terms compose the total power at the mechanical frequency /^sB(o)m) = ^^oyj^los +^sin' 
where Acos = |^+|cos((p+) + |A_| cos((p_) and Asm = |A+| sin((p+) - |A_| sin((p_). 

Given a mechanical system which is oscillating coherently at a frequency co^ (Po is simply a complex number) the single 
sided spectral density of the power at the detector, as a function of the laser detuning A, and frequency CO, will be given by 



^Pp(co,A) = n^cogK^|ain,o|'^ X 



(/A + k/2)(/(A-co^) + k/2) 



x5(co-co^) 



(A13) 



in,0 



(A2 + (k/2)2)((A-co^)2 + (k/2)2) 



x5(co-cO;^). 



For mechanical systems undergoing random oscillations, the important quantity is the power spectral density of the detected 
signal. Since this is calculated from the autocorrelation functions, it will contain products of the form Pq(^)Po(^'). Classically, 
these averages may be calculated from the Boltzmann distribution, and can be replaced with ht = ksTij/hQdm, with ks the 
Boltzmann constant, and Tij the bath temperature. Additionally, since the measured sideband is blue of the pump frequency 
(coi + co^), from the quantum theory |[23l|24]|, and the derivation below, the proper ordering to be used is the normal one (b^b), 
and thus the expectation values may be replaced with n, the number of phonons occupying the mechanical mode. As such, we 
can effectively use the derivations shown above, in both the classical and quantum cases, substituting Po by V^, and replacing 
the delta functions 5(co — co^) with unit-area Lorentzian functions. A fully quantum mechanical derivation of this result is shown 
below. 



2. Quantum Mechanical Derivation of Observed Spectra 

In this section we use the following conventions for Fourier transforms and spectral densities. Given an operator A, we take 

A{t) = ^ Jco^-^^^A(co), 
i(0)) = e''''A{t), 

J CX3 

Additionally we define the symmetrized spectral density as 5^a4(o)) = ^{Saa{(^) +^A4(— co)), and one-sided spectral densities 
which are those measured by the spectrum analyzer as 5a(co) = 2Saa{(^)- Starting from the quantum-optical Langevin equations 
for the mechanical (b) and optical (d) annihilation operators, 

b{t) = - (^io^m^^^b -igd'^d - ^/Yibin and (A14) 
d{t) = - (^/A+ I) <^ - igd{P^b) - y/Ke/2din{t) - V}ddin^i{t), (A15) 
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we linearize the equations about a large optical field intensity by displacing a 
are then given by 



' ao Then in Fourier domain the fluctuations 



S(co) 
a(co) 



-Vl^^in(co) /G(a(co) +a'^(co)) 



/(co^-co)+y,-/2 

-VK^^in(CO)- 



/(co^-co)+y,-/2 ' 



/(A-co)+k/2 



(A16) 
(A17) 



where k' = K — Kg/2 denotes all the optical loss channels which go undetected (i.e. information is lost) and G = gao. For an 
ideal measurement, Kg/2 = k, and k' = 0, so the intrinsic vacuum fluctuations (ain,0 never enter the optical cavity. For a double 
sided coupling scheme, such as the one with a fiber taper, K = K^ + K^, and so k' = Kg/2 at best, due to the back reflection from 
the cavity, which contains information about the mechanics which is lost. 

We account for all the fluctuations (vacuum and thermal) incident on the photodetector, and calculating the spectra of each 
term, we find the heterodyne detected signal. Using Equations ( A16| ) and (A17) we arrive at the operator for the mechanical 
fluctuations. 



/W^in(CO) 



/(co^-co)+y/2 



iG v'Kg/2ain(co) + VK^ain,i(co) 
/(A - co) + k/2 /(co^ - co) + y/2 



(A18) 



-/(A + co) + k/2 /(co^ - co) + y/2 
where 0)^ is now the optical- spring shifted mechanical frequency, and y = y +yoM, the optically damped mechanical loss-rate. 



a. Simplified Result under RWA ( K^/ 16co^ ^ 1 j and Weak-Coupling f G <C Kj 
Assuming that A = 0)^ and that we care mainly about the system response around co^ (where ^(co) is peaked), the relation 



( A18 ) can be simplified. 



^(co) = - v1^^.n(co) ^ 2iG 7W2^. + V^.- (CO) ^ / G 
^ ^ /(co^-co)+y/2 K /(co^-co)+y/2 V^co^ 



and we drop the term oc (RWA). 

We find using the input-output boundary condition 



^out(CO) = ain(CO) + VK^a(C0)+£:Lo8(C0) (A20) 
= ain(CO) 1 



K K 2k /(co^ - co) + y/2 



2K'Kg 4|GP /K'Kg 1 



K^ K V 2k2 /(co^ - co) +y/2 



+Sin(CO) (/Ga/^- ) +£lo8(C0) 

' ' \^ V k2 i(co^-co)+y/2 J 

= ^ii(C0)ain(C0) +^opt(CO)^in,K^) +'^12(C0)Sin(C0) +£lo8(C0) (A21) 

with the scattering matrix elements above defined as in Ref. |9 |. For the case where the mechanical bath is at zero temperature, 
the spectral density will be given simply by |^ii(co)p + |^opt(ti^)P + |'5'i2(o))P = 1, as a result of all input fluctuations being 
uncorrected, and therefore no feature is present at the mechanical frequency. 
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For the case of the n^, > 0, we find the autocorrelation of the detected normalized photocurrent /(f) = flout (0 +<^out(0 to 
5//(co) = |*ii(co)|2 + |nopt(co)p + |si2(co)p(ni, + l) + |^i2(-co)pn;, 

= l+nfe(|*i2(C0)|2 + |*i2(-C0)P) 

^2k K l^(0)™-C0)2 + (Y/2)2^(C0^ + C0)2 + (Y/2)2 

where for the last step we've used the fact that in the highly sideband-resolved regime, n = JiUij/y. 



3. Optomechanical Damping 



For the case where A = co^ (pumping on the red side of the cavity), the steady- state phonon amplitude can be written for a 
sideband resolved system far from strong coupling as (231 EH 

nb, (A23) 



yoM+Y/ 



with niy the equilibrium mechanical mode occupation number determined by the mechanical bath temperature, the mechanical 
coupling rate to the bath and 



YoM , (A24) 



the resonant optomechanical damping rate. 



Appendix B: Fabrication 



The nano-beam cavities were fabricated using a Silicon-On-Insulator wafer from SOITEC (p = 4-20 ^1-cm, device layer 
thickness t = 220 nm, buried-oxide layer thickness 2 ^m). The cavity geometry is defined by electron beam lithography followed 
by inductively-coupled-plasma reactive ion etching (ICP-RIE) to transfer the pattern through the 220 nm silicon device layer. 
The cavities were then undercut using a HF:H20 solution to remove the buried oxide layer, and cleaned using a piranha/HF 
cycle 1.31 J . The dimensions and design of the nanobeam will be discussed in detail elsewhere. 



Appendix C: Device Parameters 

Under vacuum and cryogenic parameters, the optical resonance was found to have 2^ = 4 x 10^ (corresponding to K = 
500 MHz), 0)^/271 = 195 THz (corresponding to Xo = 1537 nm), and a resonant transmission contrast of AT = 25%. The me- 
chanical mode was found to have = 1.06 x 10^ (corresponding to y = 35 kHz) and co^/2;r = 3.68 GHz. The optomechanical 
coupling rate was found to be g/27r = 910 kHz. 



Appendix D: Experimental Setup 

The detailed experimental setup used to measure the cooling spectra and the electromagnetically-induced transparency (FIT) 
window of the optomechanical crystal is shown in Figure ([5]). The setup is designed to measure both the FIT-like reflected signal 
and the transmission signal of the laser used to cool the mechanical system (though not simultaneously). 

As a light source we use a fiber-coupled, tunable, near-infrared laser, (New Focus Velocity, model TLB -6328) spanning 
approximately 60 nm centered around 1550 nm, which has its intensity controlled by a variable optical attenuator (VOA). 
A small percentage (10%) of the laser intensity is sent to wavemeter (WM, High Finesse, WS/6 High Precision) for passive 
frequency stabilization of the laser. To minimize polarization dependent losses on the electro-optical-modulator (EOM), a fiber 
polarization controller (FPC) is placed before it. 

The EOM is driven by the microwave source (RF S.G., Agilent, model E8257D-520). The RF signal is composed of an 
amplitude modulated RF- signal carrier swept between A = 1 — 8 GHz modulated at the lock-in detection frequency, cOli. As 
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Cryostat 



FIG. 5: Expanded experimental setup to include optical switches SWl, SW2 and SW3. The blue lines indicate the optical path for the cooling 
measurement (the '0' position of each of the switches), while the dashed black lines indicate the alternative switched paths (the '1' position 
of each of the switches). A single tunable laser is used as the cooling laser and mechanical transduction laser. A wavemeter (WM) is used to 
track and lock the laser's frequency to an absolute and relative value better than 100 MHz and 5 MHz, respectively. A calibrated (to better 
than 0.01 dB) variable optical attenuator (VOA) is used to set the cooling laser power. The transmitted component of the cooling laser beam 
that is sent into the optomechanical cavity is directed to an erbium doped fiber amplifier (EDFA), where the optical signal is pre-amplified 
before being detected on a high-speed photodetector (D2). The measured photocurrent from D2 is sent to a real-time spectrum analyzer, where 
the mechanical noise power spectrum is measured. A slowly modulated optical signal, on-resonance with the optical nanobeam cavity, is 
generated from the cooling laser beam via an amplitude electro-optic modulator (EOM) driven by a microwave source (RFSG). The reflected 
component of the on-resonance laser signal injected into the cavity is separated from the input via an optical circulator, sent to a photodetector 
(Dl), and then demodulated on a lock-in amplifier. Paddle-wheel fiber polarization controllers (FPCs) are used to set the laser polarization at 
the input to the EOM and the input to the optomechanical cavity. 



a result, the EOM modulation produces two probe sidebands at ±A, each with a small amplitude modulation at the lock-in 
frequency. 

A small portion of the signal from the EOM output (10%) is used as a DC control signal to keep the EOM level locked, 
compensating for any low frequency power drift during the experiment. The remaining laser light is passed through a circulator, 
a switch (SWl), and a fiber polarization controller (FPC). It is then coupled to a tapered and dimpled optical fiber (Taper) which 
has its position controlled with nanometer- scale precision. 

Switches 2 and 3 (SW2 and SW3) determine the path that the light transmitted through the taper follows. In the normal 
configuration, the transmitted light is optically amplified by an erbium doped fiber amplifier (EDFA) and then detected by a 
high-speed photoreceiver (D2, New Focus, model 1554-B) which is connected to a real-time spectrum Analyzer (RSA, Tektronix 
RSA3408B). Detector 3 (D3) is used to measured the DC transmission response of the cavity. All the other configurations are 
used to calibrate the system as discussed in the calibration section. 

Any reflected signal coming from the taper/device is detected by a high-gain photodetector (Dl, New Focus, model 1811) and 
its signal is sent to a lock-in amplifier (L.I., SRS-830). The output from the in-phase and quadrature signals from the L.I. are 
recorded, producing the reflection scan shown in Figure (3c) of the main text. 



Appendix E: Calibration for Mechanical Mode Thermometry 

To perform accurate mode thermometry, several optical switches were incorporated into the setup as shown in Figure ([5]). A 
2x2 switch (SWl) was positioned on the input/output ports of the fiber taper to control the direction of light through the taper, 
allowing the characterization of taper insertion loss asymmetry. Another 2x2 switch (SW2) was placed at the input/output ports 
of the EDFA, allowing the characterization of the optical gain, Gedfa- Lastly, a 1 x 2 switch (SW3) was inserted before the RSA 
to switch the optical path between the RSA and a power meter (D3) that reads the total power, Prsa, reaching the RSA, which 
allow us to monitor total insertion loss and provide a calibration for the electronic gain, G^. 

At the beginning of a measurement, the power into the taper from both SWl paths is measured (Pq and Pi for the '0' and '1' 
position of the switch respectively). The insertion loss of the taper Ltaper is also measured. Finally, a pickoff power is recorded 
for each of the SWl paths providing a correspondence between the measured taper input powers and the RSA powers, Prsa,o 
and /^sA,i • From these calibration values, the insertion loss before (Lq) and after (Li) the device (when the fiber taper is coupled) 
can be computed by assuming the bistability shift of the optical mode when sweeping the laser from blue to red is proportional 
to the dropped power. If we let primed values represent measurements made with the taper coupled to the device, and let AA- be 
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the bistability shift, we have 



Pin,0 PqLq AXq 
Pin,l PiLi AXi 

p' p' 

J J _ J ^RSA,0 _ J ^RSA,1 

^0^1 — Maper ~ — Maper ~ (i^^ j 

^SA,0 ^SA,1 

from which Lq and L\ can be extracted. As a consequence, the intracavity power, Pin, can be accurately determined. 

We directly measure Ge, by setting the attenuator to dB attenuation to measure Prsa and the corresponding DC bias voltage 
Vdc, and defining Ge = Vdc/^rsa- The purpose of this is a technical one; the dynamic range of an optical power meter is much 
larger than that of a voltmeter. This allows us to accurately determine the DC bias voltage of the detector for any amplitude of 
optical signal through 



This value is critical because for photodetector signals of the form VpD{t) = V{^q{1 -\-^smQ.t) where P is the modulation depth 
and CI is the modulation frequency, we have simply 

P=^^, (E4) 

^DC 

where Pq is the integrated spectral power at Q,, and we have assumed a detector load of Rl and an RSA that reports Vrms- 

During the measurement, two calibration values are measured for each point in the power-dependent cooling run: Ggj^p^ and 
^RSA- EDFA gain is measured by utilizing SW2 to insert and remove the EDFA from the optical train, while measuring 
a fixed tone at the mechanical frequency co^ generated by the Electro-Optic Amplitude Modulator (EOM). The ratio of the 
integrated spectral power of the tones gives Ggj^p^. We also measure / rsa without the EDFA in line so that when the EDFA is 
included, we have instead of Equation ( [E4| ), 

V2PqRl 

P - T^f — 

^EDFA^e^RSA 

to account for the additional optical gain. 

Finally, by integrating the Lorentzian component of the power spectral density from Equation ( A13 ), we relate the detected 
signal on the spectrum analyzer to the calibrated values shown in Equation ( [E5] ). This allows us to make accurate determinations 
of the phonon number and mode temperature. 



Appendix F: EIT Measurements 



Here we will show how the amplitude modulation of the signal sideband A is used to measure the reflection (|r((jL))p) of the 
signal reflected from the cavity. The output of the EOM can be written as: 



<^out(0 = ^in 1 + P(^l +mLlCOs(COLlO) COs(A^) 



(El) 



where the input field amplitude ain{t) = ao cos{(£>it) , ao = ^/P^Jhx^\, P is the EOM-modulation index, cOli and mu are, re- 
spectively, the frequency and amplitude modulation index on the RF signal at A. For the measurements shown in the main text 
mLi = 1. In this case one can write the field of the EOM output (cavity input) in the time domain as: 



^out(0 ~ 



cos(coi^) + ^ cos((coi + A)r) +cos((coi - A)r) 



cos((coi + A + coli)0 + cos( (coi + A - coli)0 
cos((coi - A + coli)0 +cos((coi - A-coli)0 



(F2) 
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The reflected signal is filtered by the cavity dispersion and considering the case where the pump is on the red- side of the cavity 
(coi < cOo) the reflected field is: 



aR{t) = r(co^) ^ cos((coi + A)t) + cos((coi + A)t + (coli^ - cp)) + cos((coi + A)t - (coli^ - (p)) 



(F3) 



First we assume that r(co^) is roughly constant over a range of cOli which is true for cOli < (Y; +yoM)/2. This implies that the 
smallest transparency window we could measure is on the order of the lock-in detection frequency. This limit on the transparency 
window size is reflected on Fig. 4c, where only transparency windows larger 200 kHz are reported. 

We can now write the time average detected power spectral density on the photoreceiver (Dl in Figure ([5])) by taking the 
absolute square value of the reflected field and keeping only the terms with frequency smaller than the detector bandwidth. In 
this case: 



3 + 4cos(coli^ - (p) + ^ cos(2coli^ - 2(p) + 0(2coi)] 



where /?pd = 1 AAV is the detector responsivity, Gpd = 40, 000 V/A is the detector gain and Rl = 50 is the load resistance. 

This signal is then sent to the lock-in which can measure independently the in-phase (X) and quadrature (Y) power spectral 
densities at COli: 

YU, = ^^^&^k(co.)Psin(cp) (F4) 
It is then easy to see the reflection amplitude and phase are given by 



a^p^7?pDGpD V Li Li X I COli 

From the imparted change in the phase the signal delay is then calculated as: 

^(R) ^ _1_ 
COli 

where x^^^ > (x^^^ < 0) represent a delay (advance) on the signal. 

Here we have neglected the gain provided by the lock-in, which is important to determine the absolute value of r(co^). To 
account for that we calibrate the X channel by a normalized transmission curve taken with low input power. Our assumption is 
that the cavity-taper coupling is not affected by the input power. A analogous result can be found for the case where the control 
laser is on the blue side of the cavity (coi > 0)^). A more detailed description of the FIT experiment can be found in Ref. 1261 . 



Appendix G: Analyzing the Mechanical Mode Spectra 

To determine the total spectral power at co^ for a given measured spectra, we first subtract a background taken with the cooling 
laser far-detuned from the cavity (in the same calibration conditions). We then perform a least squares fit to a Lorentzian function 
of the form 



L(co) = -2 (Gl) 

(^) +1 



with fit parameters A, co^ and y. The spectral power is then given simply by 

^co„ = ^. (G2) 

To extract the intrinsic linewidth we first fix the input power Pin. We then lock the pump on the red side of the cavity (at 
A = +03^) and measure the total linewidth, yred = 7? +7om^- We repeat the measurement on the blue side (at A = — co^), where 
ybiue = Ji ~7oM^^- Using low input powers where y^^^^^ ^ 7' avoid amplification of the mechanical oscillations, we have 
TSM'=7gM^ which leads to 

y.= M±y^, (G3) 
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Appendix H: Measurement of Phonon Number and Error Analysis 



Equation ( A13| ) shows an explicit form for the sideband power ampHtude seen by a photodetector for a red detuned pump 
laser. More specifically we can find a relation between the number of phonons inside the cavity and the power spectrum for our 
experimental setup. As shown before, the RF-spectra are detected via a RSA which displays the power spectral density of the 
voltage comi ng fro m the photodetector (with gain Ge and Gedfa)- The single sided power spectral density at the detector is given 



by Equation ( A13 ), and denoted 5'pp(co). Since the power is related to voltage by an electronic gain G^, then Svv{(^) = G^5'pp(co). 



When the EDFA is used, there is an additional gain factor, and 5'yy(co) = Ggj)p^G^5'pp(o)). Finally, the RSA reports power as 
opposed to squared voltages, and so the final spectral density measured is 5'(co) = 5'yv(co)/2/?z,, where Rl = 50 Q, is the input 
impedance of the RSA and the factor of two in the denominator comes from the conversion of peak-to-peak voltage to RMS 
voltage. Then, in terms of integrated power, the power detected in the sideband on the RSA, Pco^ is related to the heterodyne 
detected integrated spectral density by the relation 



(Gg Gedfa )^ 
2Rl 



(HI) 



We would like to write this equation as a function of all the independent variables measured for the system, and from that 
estimate the error on the measured number of phonons. 

From the DC transmission spectra, the optical components K, and cOq can be determined. Both Gedfa and Ge are measured 
and latter compensates for any discrepancy in the value of Rl. From the RF-spectra one can determine the total mechanical 
linewidth y = y +yoM, the mechanical frequency co^ and the total RF-power Prsa- The FIT spectra give the true detuning A 
between the pump laser and the cavity and, as shown before, can be used to determine the power at the cavity Pi^. 

Using Equations ( |A24| ) and ( |A9| ) we can rewrite the integrated form of Equation ( |A13| ) as: 



(K./2)Mr 



(A-co^ 



-{K/2y 



K(y-y-)/2 



(H2) 



From Equation (|H2|) we can write the expression that relates the number of phonons, n, and all the system parameters as: 



f 2Rl 



1 



(A-C0^)^ + (K/2)^ 
(K,/2)Pin 



(H3) 



We can calculate the cumulative error for the number of measured phonons on the cavity using this equation and based upon 
on the measurable variables which gives: 



An 
n 



(y_y.)2 y_y.)2 



(H4) 



k/2 



(K/2)2 + (A-C0^)2 K 



5k^ 



2(A-co^) 



(K/2)2 + (A-co^)2 



5A^ 



2(A-co^ 



(K/2)2 + (A-co^)2 



1/2 



Here we neglected the error on Ge and Gedfa which are much smaller than any other error quantity. To determine the variation 
for K, Ke and COq, we measured the DC optical spectrum for every single data point in Figure (4a) of the main text and determined 
5k, dKe and 5cOo from the normalized standard deviations of each of the values. The measurement uncertainty of these values 
are below 0.7%. The mechanical properties 5y, 5Prsa and 5co^, were determined from the deviation on the spectra fits using 
a 95% confidence interval, which produces percent errors below 0.6%. The pump laser detuning from the cavity is controlled 
by the EIT reflection spectra. To find the variation of the detuning 5A we once again computed the standard deviation of all the 
measured detunings, which results in a deviation of less than 0.3%. 

Finally, the two main sources of error in our data are the determination of the intrinsic mechanical quality factor (reflected in 
y) and the input power, Pin. The uncertainty in the mechanical linewidth, 5y , is found by repeatedly measuring it at a single 
power level and computing its standard deviation (found to be ^ 1.6%). Using the calibration procedure discussed above for 
Pin, the error lies in the determination of losses Lq and Li . In the worst case the calibration would be off by the ratio between the 
input loss, Lo in the present experiment, and the square root of total loss y/LoLi producing a percentage error of ~ 4.0% to the 
input power. 

Taking all of these factors into account produces an overall uncertainty of ~ 4.5% in the measured absolute phonon number. 
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Appendix I: Estimating the Temperature Shift from Thermo-optic Effects 

Absorption in the dielectric cavity causes the temperature of the dielectric cavity to increase locally. This effect is expressed 
through shifts in the refractive index of the structure, and the thermo-optic coefficient of Silicon ll32ll . As such, we can estimate 
the temperature of the cavity by looking at the shift in the cavity frequency, starting from a known temperature. 

The starting point of the analysis is the cavity -perturbation formula for dielectric cavities [33J, 



co-coq l /6£(r)|E(r)pdr 
^o^^"2 /£(r)|E(r)pdr • 



(11) 



From the relation e/eo = we find 58 = 2ndmo. By assuming that the cavity as a whole is heated to a temperature Tq, the 
integral in Equation ([iT]) can be written as 



co-o)o ^ -^(ro)coo- 



/si|E(r)pdr 



7K7b))2|E(r)Pdr 
Using the values of ^(r) found in literature 1321 , and a value of 

/si|E(r)pdr 



x{n{T)-n{To)). 



(12) 



/K7b))2|E(r)Pdr 



7.5066 X 10- 



calculated from the finite element simulations (FEM) of the mode profiles, we plot the wavelength shift from 17.6 K up to 300 K 
in Figure ([6^). The total shift of 12.5 nm agrees with the experimentally observed change in resonance wavelength. 
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FIG. 6: a, the measured wavelength shift compared to the theoretical shift predicted by Equation (|I2]) for a range of cavity temperatures, using 
17.6 K as the reference point, b, the measured power-dependent wavelength shift of the cavity with fitted individual contributions due to free 
carrier dispersion (blue) and refractive index change (red), as well as their sum (black), for the two bounds discussed in the text. 

This analysis can be applied to the wavelength shift data for various input powers at low temperature where the initial cavity 
temperature is measured by thermometry methods discussed above. We note an initial blue-shift of the cavity, which is attributed 
to free-carrier dispersion effects 1 34] and can be modeled by a power law dependence on intracavity photon number, A{nc)^. The 
temperature-dependent data for the refractive index of Silicon in Ref. 1 32 1 is valid only for T > 30 K so the power-dependent 
cavity heating for a starting temperature of 17.6 K, for the largest intracavity photon number, can only be bounded. For the 
upper bound, we assume dn/dT = for T < 30 K, resulting in a AT^iax of 16.8 K. For the lower bound, we assume dn/dT = 
dn/dT\T =30 K for 7" < 30 K, resulting in a Armin of 7.8 K. These bounds and their respective fits are shown in Figure 



Appendix J: Temperature-Dependent modifications to the intrinsic mechanical damping 

Independent measurements of the mechanical quality factor, Qm, at varying bath temperatures indicate that the Qm changes 
with temperature (Figure (|7^)). These measurements are taken at low intracavity photon number, rendering free-carrier effects 

negligible. As such we can model the mechanical loss rate as JiiT) = Ji^riT) -^J^ where y|^^ is the measured loss rate at the 
reference temperature (17.6 K). The extracted form of Ji^riT) is shown in Figure (ItI)). 
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FIG. 7: a, the measured intrinsic mechanical quahty factor for various cavity temperatures, b, the inferred intrinsic mechanical loss rate due 
to temperature, y/ modeled by a polynomial fit. 



Appendix K: Photon-number dependent modifications to the intrinsic mechanical damping 

The deviation of the expected cooled phonon number from the measured value is a result of two factors: bath heating and 
an increase in the intrinsic mechanical loss rate (y^) due to heating and free carriers. Since the integrated spectral power of the 
mechanical mode depends only on the product y^T^ (for large intracavity photon numbers ric), naively ignoring the latter effect 
results in an estimated change of AT > 50 K in the bath temperature for 2000 intracavity photons. This is unrealistic as such a 
temperature change would tune the optical mode red by > 300 pm (from Equation (|I2])), while the actual measured shift is closer 
to 10 — 20 pm (from Figure ^)). In fact through independent measurements (where dynamic back-action was minimized), we 
found that the mechanical linewidth is a function of the number of photons in the cavity. We attribute this to a nonlinear process 
in the cavity involving the generation of free carriers, which will be explored in depth elsewhere, and introduce an additional 
loss channel y^^pc in the mechanical loss rate so that we have 

y.- ^ Ji = yf^ +y-,r(rK)) +y-,FcK). (ki) 



From the relations shown on previous sections, we have J^Jqi^^ = y|^^ +yoM. Incorporating Equation dKll), we have experi- 



mentally, ycooied = y +yoM, with their difference yielding the magnitude of the additional loss rates. However, for A = co^ and 
He > 10, yoM tends to be large compared to y^, making this subtraction quite error prone. Thus, to get accurate data for high 
intracavity photon numbers we use a range of larger detunings, noting that yoM for A co^ and fixed ric (approximately). 

This loss is then modeled using a power law dependence on ric (Figure ([8^)). 

Using the models of mechanical loss from Figure ([8^) and Figure ([T])) with the thermometry technique outlined earlier (making 
the replacement to yO allows a more accurate determination of the temperature rise in the cavity, as well as the characterization 
of the individual contributions of Jij and y^^pc as a function of ric. The result is an estimated increase of 13.2 K in at the 
highest input power, well within the previously fitted temperature bounds. 

The addition of a free carrier related loss channel is further corroborated by pumping the Si sample above the band gap with 
a 532 nm solid state green laser, directly stimulating the production of free carriers. The degradation in can be only partially 
explained by heating due to absorption since the maximum 19 K temperature rise estimated from the cavity red-shift results in an 
expected of approximately 70,000 at the highest power (Figure (|7^)), whereas a far lower value is measured. The remaining 
excess loss is attributed to the presence of free-carriers. These results are shown in Figure ([9]). 



Appendix L: Shot Noise Considerations 

We consider here the impact of using a non-ideal amplifier (EDFA) on the measured signal, and the deviation from quantum 
limits. For a coherent optical beam with frequency coi and power P incident on a photo detector, the single sided power spectral 
density of the shot noise is simply 

^shot(co) = ^/2m^, (LI) 
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FIG. 8: a, excess loss as a function of ric, inferred from far-red-detuned (A > 5.5 GHz) measurements of Ycooied- b, breakdown showing the 
individual contributions of y-^^ (gray), jij (red), and y/^pc (blue) to the total y/ (o). 
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FIG. 9: The Qm degradation as a function of 532 nm laser power. The purple line shows the expected Qm for the bath temperature rise inferred 
from the wavelength shift data. The deviation of the data from this prediction suggests an additional loss channel related to the presence of 
free carriers. 



independent of frequency. As part of the calibration procedure, Gedfa is measured to characterize the gain provided by the 
EDFA and Ge is measured to characterize the transimpedance gain and quantum efficiency of the photodetector. These values 



can be used in conjunction with Equation (LI ) to predict the expected spectral background assuming only the presence of shot 
noise and noise-free gain. To wit, 



^tr'^'^ = A/2to.G2 (L2) 



where the prime indicates the value has been adjusted to account for insertion loss from the cavity to the photodetector. The 
difference between this predicted level and the measured background level, 5'background. gives the non-ideality of the EDFA and is 
attributed to an excess noise which also includes the amplified spontaneous emission (ASE) noise 05 i . This deviation is shown 



in Figure ([TOp). Defining 

- s - . . - 

-'background V'^shot 



r.2 — ^2 _ /^(amplified) N 2 /j ox 

"^excess "^background Wshot / v^-^/ 



this additional noise reduces the measured signal-to-noise ratio (SNR). We can predict the shot noise limited SNR, representing 
the largest measurable SNR for the current experimental setup set assuming perfectly efficient detection from the measured 
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FIG. 10: a, a comparison of the measured background spectrum against the shot noise level amplified by an ideal, noise-free amplifier, b, a 
comparison of the measured signal-to noise ratio (SNR) to the maximum achievable SNR for the experimental setup assuming an ideal, noise- 
free amplifier and perfect quantum detection; the predicted SNR uses only measured device and calibration parameters, along with ^'g^cess 
determined from (a); the purple dashed line shows SNR assuming an overcoupled system, where Kg = K, with noise-free gain and the black 
line shows the same system without taper/insertion loss (representing the most ideal of ideal cases). 



device and calibration parameters using Equation (H2 ) and ( |L1| ). This ratio is given by 



2to,/^; • ^^^^ 

Similarly, the expected measurement SNR, using only device and calibration parameters, is given by 

_ 4GiDFA^^-PsB/(y^+yoM) 

^JNKpredicted - 2 r2pf 1 c2 ' ^^^^ 

^''^^o^EDFA^^eMn ^ ^excess 



which closely corresponds to the measured values, as seen in Figure ( 10 3). 



